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In this paper, two kinds of analytic singular solutions (finite-time and infinite-time singular 
solutions) of two classical wave equations (the Boussinesq equation and a generalized KdV equation) 
are obtained by means of the improved homogeneous balance method (HB method) and a nonlinear 
transformation. The solutions show that special singular wave patterns exist in the classical models 
of shallow water wave problem. 
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(NJ . I. INTRODUCTION 

All kinds of wave patterns [1], from linear systems to nonlinear systems, contribute to understanding complex wave 
phenomena. Generally, the wave patterns are symmetric resulting from the fascinating effects of instability. When 
the instabilities of a nonlinear system are violent, the energy may focus into a spike at a point where special wave 
patterns (singularity in the nonlinear system) appear; the special wave patterns are local divergences in the amplitude 
or gradients of some physical field. As is known, a lot of spatio-temporal systems exhibit this phenomena, e.g., in the 
turbulence problem it is observed in Ref [2] that large fluctuations in the derivatives occur and the dissipation field 
f""^ | appears multifractal (a set of nested singularities); in the surface wave problem, the singularity formation has been 
observed and studied [3,4]; nonlinear optical systems exhibit self-focusing effects, which may lead to the collapse of 
the optical power density into local divergences which may have important consequences on the integrity of optical 
fibers and laser systems; etc. 

For understanding singularities of nonlinear physical systems, it is helpful to study singular solutions of partial 

1 " i l differential equations (PDEs) that model nonlinear physical systems, and many papers have been written on the 

problem; for example, near-singular solutions of the Navier-Stokes equations and singular solutions of the Euler 
equations [5]; the possibility of singularities arising in Burger's equation [6]; some methods were presented for the 
investigation of the singularity formation in the NLS equation [7] ; near-singular solutions of the complex Ginzberg- 
Landau equation from the view point that the NLS equation is the conservation limit of the complex Ginzberg-Landau 
equations [8]; etc. However, most methods used in studying singular solutions are based on perturbation techniques 
' or the singular point analysis. Also, the explicit forms of singular solution are seldom discussed, except for rational 
solutions. In this paper, we try to obtain the analytic singular solutions of nonlinear PDEs by two direct methods: an 
t-H ■ improved HB method [9] and the invariant-Backlund transformation [10] based on a special nonlinear transformation. 

The homogeneous balance method [11] has shown its efficiency in finding analytic solitary wave solutions of many 
£H . PDEs. Its essence can be presented as follows: the nonlinear PDE is given by 
• i— I ' 

£H \ P{u,U x ,U tl U xx ,U xt ,UUi ) = 0. (1) 

Supposing the solution of equation (1) is of the form 

N 

u = £>/ W *)). (2) 

where i is an integer, a; are constant coefficients, f(tp(x,t)) is a function of (f(x,t), ip(x,t) is a function of x and 
t, and superscript (i) represents the derivative index. According to the assumption of homogeneous balance, in 
the PDE (1), the nonlinear terms and the highest order partial derivative terms ought to be partially balanced. 
Then N is obtained, the expression of f{tp{x,t)) and the relation of /W . fU) and ) can be derived. Assuming 
ip(x, t) = 61n(l + e ax+,3t+1 ) and substituting formula (2) into equation (1), after deciding coefficients a%, b, a, (3, 7, the 
solitary ware solution of the nonlinear PDE is obtained. Recently an improved HB method [9] has been reported: 
after balancing the nonlinear and the highest order partial derivative terms, substituting f(ip(x,t)) and the relation 
of jW . fU) and j n t equation (1), we get 

F(f, /", ...ip x , tp xx , . . . ipt,(p tt , ■ ■ ■ Vxt, ■ ■ •) = 0, 

where F is a function of /',/",.•• <Px> <Pxx, ■ ■ ■ tfu Ptt> ■ ■ ■ <Pxt> ■ ■ ■ Obviously, F is a linear polynomial of /',/",.■ ■ Setting 
the coefficients of /',/",.. . to zero yields a set of partial differential equations of (p(x, t). Taking note of the terms of 
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/', their forms must be f' <p Xl , X2 ,... Xp , tl ...t q - So the coefficient of/' in F(f',f", . ..<p x ,tp xx , . ..<p t ,<Ptt, ■ ■■<Pxt, ■ ■ •) = is 
k 

a linear polynomial a i<P X i, X 2,...x pi ,ti...t<,i, an d then the set of partial differential equations has an important feature, 

i.e. the equation from the coefficients of /' is a linear PDE of ip(x,t), whilst the other equations can be regarded as 
constraint conditions to the linear equation. Thus solving the set of PDEs reduces to solving the linear PDE with 
some constraint conditions, and then the analytic solution (including travelling and non-travelling wave solutions) of 
nonlinear PDEs can be obtained by the improved HB method. 

As an example, the Boussinesq equation and the KdV equation (fundamental models of the shallow water wave 
problem) are discussed. Analytic singular solutions of the equations are obtained by the improved homogeneous 
balance (HB) method. An particular, a finite-time singular solution of the Boussinesq equation is obtained in this 
paper, which was produced from a non-singular physical field in the process of time evolution. The finite-time singular 
solution is relevant to a conjecture — in hydrodynamics many authors have shown that singularities of fluid motion 
can be formed in finite time, and some simplified models of fluid motion have been presented to demonstrate the 
conjecture [12]. Here the Boussinesq equation gives another simplified model of fluid flow to do the same job, since 
this equation is a fundamental model of the shallow water wave problem. 

The Backlund transformations were developed in the 1880s for use in the related theories of differential geometry and 
differential equations. Afterward, the relationships between the Backlund transformations and the inverse scattering 
transform [13] (1ST) or the bilinear form [14] were presented, and the Backlund transformations have been used to 
find analytic solutions of nonlinear PDEs. In this paper, based on a special nonlinear transformation, the invariant- 
Backlund transformation of a generalized KdV equation is obtained and the analytic singular solution of it is obtained. 

This paper is organized as follows. In Sec II, an analytic singular solution of the Boussinesq equation is obtained 
by the improved HB method, which shows that the finite-time and the infinite-time singularity both exist in the 
Boussinesq equation. In Sec III, an analytic singular solution of the KdV equation is obtained by the improved 
HB method and an analytic periodic singular solution in space of the generalized KdV equation is obtained by the 
invariant-Backlund transformation. In Sec IV, the main results are given. 



II. BOUSSINESQ EQUATION 

The Boussinesq equation [15] is a classical model of long wavelength hydrodynamic waves and other physical systems 
and is written in the form 

utt - u xx - a(u 2 ) xx + bu xxxx = 0, (3) 
where a, b are real constants. It can be derived from the incompressible fluid equations 

d t u = -^d x p', 

d t w = -\d zV ', (4) 
d x u + d z w = 0, 

where u, w are velocities in the x, z directions, p' is the relative pressure. Eq. (3) can be derived from the Toda 
Lattice also. The Boussinesq equation (for wave propagation in both directions) describes wave motion in weakly 
nonlinear and dispersive media, and a suitable approximation enables the KdV equation (for which wave propagation 
is restricted to one direction) to be derived from this equation. Much work has been reported for this equation, for 
example, the soliton and the Painleve expansion [16], the soliton and the bilinear form (Hirota's method) [17], periodic 
wave solutions and the structure of the rational solution to Boussinesq equation [18]. . ., etc. 

Here, the improved HB method is employed to obtain the singular solution of the Boussinesq equation. Supposing 
the solution is of the form 

N 

« = (5) 

where N is an integer, substituting formula (5) into equation (3), and balancing the nonlinear term a(u ) xx and the 
linear term bu xxxx , we get N — 2, and then 

u = f ul + f uj xx . (6) 

Substituting formula (6) into equation (3), and collecting all homogeneous terms in partial derivatives of tv(x, t), we 
have 
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[fe/W - 2a/'"/" - 2a/"/ (4 Vx + [W (5) - 24a/"/'" - 2a/'/ 4) R^+ 
[/ (4) ^ - /< 4) $ + (456/W - 24a/"/" - 12a/'/>^L + (206/W- 
8a/"/" - 4a//'" V^w.^] + [(w M w 2 + w t w x w xt + w 2 w xx - 6w 2 w xa: )/" + 
(156/ - 6a/ / + (606/ - 20a/'/ )u x u xx uj xxx + (156/ - 2a/'/ ) 

^x^xxxx] ~t~ 

[(2u 2 t + 2uj x lu xU + lott^xx - 3w 2 x + 2u> t u) xx t - ^xUJ X xx)f + (106/ - 2a/ / )u xxx + 
(156/ — 2a/ / )lo xx uj xxxx + 66/ uj x uj xxxxx } + (oj xx tt + buj xxxxxx — w xxxx )f = 0. 



(7) 



Setting the coefficient of lu x in (7) to zero yields an ordinary differential equation for /, namely 

bf^ - 2a/"/" - 2a/"/ 4 ) - 0. (8) 
The solution of equation (8) is obtained as 



/ = --lnc, (9) 



a 



which yields 



f" f "' - A f (5) f' f (4) = 36 f (5) f" f" _ 6 W4) 



<■"' _ 26 i-(4) f" _ 36 J-'" f ' f ' _ 66 J-" 

J J a •' 1 •' •' a * 1 ^ a •' ' 



>' ' ,7 66 5" ' ( 10 ) 



Substituting formulae (10) into equation (7), it can be simplified to a linear polynomial of /', /", . . .; then setting the 
coefficients of /', /", ... to zero yields a set of partial differential equations for uj(x,t), 

<A - u\ - 3bu)l x + 4bu) x u) xxx = 0, (11) 
uJtt^l + 4uj t uj x uj xt + uj 2 uj xx - 6ujIuj xx - 3bu> xx + %uj 2 x uj xxxx = 0, (12) 
2 w xt + Z^xUixtt + i^tt^xx — 3lu xx + 2iu t Lo xx t — 4:LO x ui xxx ~ 2ho xxx + 2>bto xx u> xxxx + Qbto x to xxxxx = 0, (13) 

^xxtt ~t~ biO xxxxxx L^xxxx — 0. (1 4) 

We note that the equation (14) from the coefficients of /' is a linear PDE for u>(x, t). Then equations (11)-(13) can be 
regarded as constraint conditions to the linear equation (14); thus solving equation (3) reduces to solving the linear 
PDE (14) with some constraint conditions (11)-(13). 
Equation (14) may be integrated once to yield 

w x tt + bu> xxxxx — lo xxx = p(t), (15) 

where p(t) is an arbitrary function of time. It's easy to know that equation (14) has the solution 

u;(x,t) = S(Z) + q(t), (16) 

where S(£) is the traveling wave solution of equation uj xtt + bu xxxxx — lo xxx = 0, and q(t) satisfies equation ^(g(t)) = 
p(t). So the solution of equation (14) is of the form 

uj(x, t)=d Q + di(x - vt) + d 2 (x - vtf + d 3 e^ x - vt) + die-^ x ~ vt) + q(t), (17) 



where (3 — y * > 0; note that solution (17) is a non-travelling wave solution of equation (14). Substituting formula 

(17) into the constraint conditions (11)-(13), a set of ordinary differential equations are obtained. Solving this set of 
equations, we find 

d 2 = 0, (18) 
6d 2 (-l + 4v 2 ) 
3 ~ 12d 4 v 2 (-l + v 2 )' 
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Thus the solutions of equations (11)- (14) are obtained as 



U {x, t) = d + d l{ x - vt) + Z l{ -} + 'lly (X - Vt) + ^- 0{X ~ Vt) - dl{ ~ 1 + V) t, (19) 

12d,4V z {— 1 + v z ) v 

where b, v, do, ciiand di are arbitrary constants. Substituting solution (19) into formulae (9) and (6), the analytic 
solution of equation (3) is obtained, namely 

w y 1, 1 i2d 4 ^(-i+^) J 



a(d + d4 e-«— ) - + + - «*» ' 

where a, &, v, do, c?iand d\ are arbitrary constants. Here a just is a parameter to decide the relative value of u(x, t), 
and when d\ — 0, the solution (20) decays to the solitary wave solution. 

Solution (20) is very complex and includes six arbitrary constants, so it needs discussion. It's easy to understand 
that the singularity of solution (20) comes form the zero value of formula (17), namely, it is decided by uj{x, t) = 0. If 
the analytic solution x = T{t) or t = X(x) of lo(x, t) = is obtainable, the main property of the solution can be given, 
but formula (17) includes the transcendental function, so it's impossible to get an explicit solution. By numerical 
analyticity, the main property of the solution is reported as follows: when 0.5 > v or v < —0.5, b < 0, d\ < 0, d^ < 0, 
the unlimited-time blow-up solution exists, and when 0.5 > v > — 0.5(v ^ 0), b > 0, d\ > 0, d± > 0, the finite-time 
blow-up solution exists. Fig.l. shows a finite-time blow-up solution. The numerical results suggest that: (1) the 
constant v decides the region of the singularity, when 0.5 > v > 0, the singularity exists in the half-line [-co, to] and 
when > v > —0.5, the singularity exists in the half- line [to, oo]; (2) the constant a just decides the relative values of 
u(x,t), and the constants b, d , d\, d± decide the values of to, which is the blow-up time; and (3) there are two kinds 
of finite-time blow-up evolution modes at v = —0.5 or!i = 0.5 and 0.5 > v > —0.5. Fig. 2. (v = —0.5) shows a time 
evolution mode of the finite-time blow-up solution. Fig. 3. (0 > v > —0.5) shows another time evolution mode of the 
finite-time blow-up solution. 



III. KDV EQUATION AND GKDV EQUATION 

As is known, the KdV equation is not only of mathematical interest but also of practical importance. It has been 
shown to describe small amplitude shallow water waves, hydromagnctic waves in a cold plasma, ion-acoustic waves, 
acoustic waves in an anharmonic crystal, and wave motions in other biological and physical systems. The standard 
KdV equation is given by 

u t + auu x + (3u xxx = 0, (21) 

where a, j3 are real constants. When the improved HB method is applied to the KdV equation, the following analytic 
solution of KdV equation is obtained 

u(x,t) = {ASdie^-^v^dl^e^-^b + (3d 1 v(-d 1 e b(x - vt ^(8d i + 

doe b(x-vt)^ + 8d 2 b + d 2 e 2b(x-vt)( 3vt _ x ^ + ulv 2 (do + d 1 (x- ivt)))) (22) 
/{a(/3d 2 1 e 2b ( x - vt '> - 4d 4 v{d 4 + e b{ - x ~ vt \do + dl (x - 3vt)))) 2 ). 

It is easy to show that solution (22) is an unlimited-time singular solution. 

Recently, some generalized KdV equations (gKdV) have been discussed widely, including a few high order KdV 
equations [19], the q-KdV equation [20], and some generalized KdV equations [21]. One generalized KdV [10] 

u t + (n+l)(n + 2)u n u x + u xxx = 0, (23) 

has been presented and a more generalized form 

u t + a(u m ) x + b(u l ) xxx = 0, (24) 
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has been introduced by Roscnau and Hyman [22]. Equation (24) has yielded the compactons solution [23] for certain 
values of m and I, which, like solitons, have the remarkable property that after colliding with other compactons, they 
reemerge in the same coherent shape. In the paper, a concrete equation of the gKdV equation (23) (in the case n = 3) 



u t + 20u d u x + u xxx = 0. 
is considered. For a discussion of equation (25), the nonlinear transformation 

g(x,t) 



u(x, t) = 



/(*,*)' 



(25) 



(26) 



is introduced. Then, substituting solution form (26) into equation (25), and making the coefficients of like powers of 
f(x,t) vanish, we obtain the following set of equations, 



r 
r 
r 
r 
r 
r 
r 
r 



■■-§gU9g 3 + y 2 x ) = 0; 

-40 ff 3 u /x = 0; 

f (6fl 3 5x + / x 2 5x+ff/x/xx) = 0; 
-m 9 2 ulf x = 0; 
20g 2 (3u o g x + gu Qx 

-f {9 ft + Wguofx + 3g x fxx + Zfxgxx + gfx 



0; 



(27) 



c) = 0; 



60gu o (u o g x 
g t + 20ulg x 



gu 0x ) = 0; 
■ 60gulu 0x - 



where uo(x,t) satisfies the equation 



u ot + 20u o u Qx + u Qx 



0, 



0, 



(28) 



which is the same as equation (25) for u(x,t); thus the set of equations (26-28) constitutes an invariant-Backhand 
transformation of the g-KdV equation (25). 

It is obvious that equation (28) has the trivial solution u (x,t) = 0, then substituting uo(x,t) = into formula 
(26), the transformation is simply 



u = 



f(x,t)i 



(29) 



and equations (27) become 



9x 



0; 



9t 

9 ft + Zgxfxx + Zfxgxx + gfx 
6g 3 9x + fhx + gfxfxx = 0; 
9g 3 + 2/ x 2 = 0. 



(30) 



The frist and last equations are easily solved (the other equations can be treated as constraint conditions), from which 
solutions of equations (30) can be obtained. Substituting the solutions into transformation (29), the analytic solution 
of the g-KdV equation(25) is obtained as 



u(x, t) 



v 3 



sec 3 (• 



-16c+9pQ+^t) 



(31) 



where c is a constant, and v is the velocity of the travelling wave. It's easy to understand that the singularity of 
solution (31) is decided by equation tani ( ~ 16c +M? : + t ' t ) ) = q. Thus the solution (31) is a space periodic singular 
solution, i.e., the g-KdV equation contains the unlimited-time singular solution, which is periodic in space. 



IV. CONCLUSION 



In conclusion, two methods to seek analytic solutions of nonlinear PDEs are introduced: an improved HB method 
and the invariant-Backlund transformation based on a special nonlinear transformation, and two kinds of analytic 
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singular solutions (finite-time and infinite-time singular solutions) of the Boussinesq equation and a g-KdV equation 
(including the KdV equation) are obtained. A finite-time singular solution of the Boussinesq equation is discussed, 
which is relevant to singularity formation in finite time in hydrodynamics. A periodic in space unlimited-time singular 
solution of the g-KdV equation is obtained, and the invariant-Backlund transformation of the equation is presented. 
In addition, the nonlinear transformation u = -j^ + u is useful in finding solutions of equation (23). 
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Fig. Caption: Quantities plotted in the figures are dimensionless. 

Fig.l. A finite-time blow-up evolution is shown of a finite-time singular solution; the values of the parameters are as 
follows: a = —0.1, b = 1, v = —0.5, do = 0.5, d\ = 1, d 2 = 1. In the figure, the max value u(x, t) and the asymptotic 
blow-up time is indicated. 

Fig. 2. A finite-time blow-up solution is shown; the values of the parameters are as follows: a = —0.1, 6=1, 
v = —0.5, d = 0.5, d\ = 1, d 2 = 1. At t = 0, u(x,t) looks like a hump(a standard solitary wave), the height of the 
hump becomes higher and higher, and at t ~ 1.0 the height of the hump tends to infinity. 

Fig. 3. A finite-time blow-up solution is shown; the values of the parameters aret as follows: a = —0.1, b = 1, 
v = —0.2, do = 0.5, di = 1, d 2 = 1. At t = 0, u(x,t) looks like two humps, one hump runs to another hump, at 
t ~ 1.0, two humps come into collision and become a single hump, after which the height of the hump tends to infinity. 
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